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The present final part of a triad of plenary talks on the results of high energy 
\Q \ gamma-ray astronomy and the origin of Cosmic Rays is primarily devoted to the 

^ \ physics interpretation. I will start with the stereoscopic method of operating 

J> ' several imaging atmospheric Cherenkov telescopes in coincidence as pioneered by 

^ ■ HEGRA, and a summary description of the results obtained. Then I will turn to 

lO ' the search for gamma rays from Supernova Remnants over the last decade and 

^ . argue that only the quantitative comparison of observations with a consistent the- 

i ory of particle acceleration can lead to a generally acceptable picture of Cosmic 

^ ' Ray origin. The elements of this theory are outlined. It is subsequently used to 

model SN 1006 and Cassiopeia A, the two sources that until now could be inves- 
tigated to the required detail. The analysis shows for the first time that emission 



^ ■ spectra and morphological detail are in agreement with the concept that these two 

^ ■ distinctly different objects are representative members of a suspected Supernova 

source population in the Galaxy. The continuing study of this population is an 
- . . essential part of the program of the major new gamma-ray instruments. 

rS 

; 1. The HEGRA Stereoscopic System 

Representing also the HEGRA experiment in this overview I want to emphasize 
first the significance of the stereoscopic method in ground-based 7-ray astronomy. 
This observation technique positions several telescopes in the Cherenkov light pool 
on the ground so that the same gamma-ray shower in the atmosphere is observed 
from different positions, in a manner analogous to the practice of a land surveyor 
(Fig. 1). The reconstruction of the 7-ray direction is achieved on a purely geo- 
metrical basis by combining the different shower images into a single focal plane 
detector ("camera"). The shower impact point on the ground is reconstructed 
similarly from the images in the spatially seperated telescopes. An immediate 
result is increased angular resolution and a higher energy resolution. At the same 
time it is possible to strongly improve background rejection, in particular against 
local muons from unrecognized showers in the distant atmosphere. The energy 
threshold for such a system is therefore lower than for each of the individual 
telescopes seperately. 
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Fig. 1. Schematic of stereoscopic 7-ray observations. The Cherenkov light of the 
shower electrons illuminates a cone whose basis on the ground has a radius of about 
120 m. Pointing telescopes inside the Cherenkov cone see the shower under different 
angles (left panel). The elongated shower images in the focal plane (with a length 
exceeding their width) point towards the primary 7-ray direction. Bringing the 
various images into a single camera allows a geometric reconstruction of the 7-ray 
direction. An analogous image extrapolation from the seperated telescopes yields 
the shower impact point on the ground (right panel). Courtesy G. Piihlhofer. 



The stereoscopic array of the HEGRA experiment on La Palma which 
was finally dismantled in late 2002 has been combining five relatively small tele- 
scopes of 3.3 m effective diameter (8.5m^ mirror area) each. Using the technique 
described above it was able to detect and thereby to confirm all the Northern 
Hemisphere sources originally detected by the (single) 10 m (75m^ mirror area) 
Whipple telescope in Arizona. In addition, the HEGRA system detected a num- 
ber of very weak sources like the first unidentified TeV source TeV J2032+4130 
in the Galaxy, the first Northern Hemisphere Supernova Remnant Cassiopeia A 
(Gas A; Fig. 2), and the first radio galaxy in the TeV range: M87 in the cen- 
ter of the Virgo cluster. Gonsistent with their fluxes of merely several percent 
of the Grab Nebula, the three new sources could not yet be confirmed by other 
instruments. 

I shall not enter into a discussion about the overall observational results of 
TeV 7-ray astronomy which the previous speakers, T.G. Weekes and T. Kifune, 
have already presented so ably for the Extragalactic and the Galactic sources. 
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My assignment here is rather to interprete these results with regard to the origin 
of Cosmic Rays (CRs). More specifically I shall ask the question, whether and 
to which extent the TeV-detection of Supernova Remnants (SNRs) in our Galaxy 
has solved this problem which has eluded physicists for more than ninety years, 
since the discovery of Victor Hess in 1912. 




Fig. 2. Cas A observed in 6 cm radio wavelengths with the VLA (R.J. Tuffs, 1986). 
The bright ring of emission is interpreted as the result of a compressed circumstellar 
wind shell [51], with strongly amplified magnetic field [53]. 



2. Test of the Supernova Remnant origin of Cosmic Rays 

Following the hypothesis of Baade and Zwicky from 1934, that Supernova explo- 
sions might be the ultimate sources of the CRs [1], the idea has been entertained 
in various forms. Indeed as time progressed it became more and more clear that 
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there are not many viable alternatives from an energetics point of view. And it has 
slowly become a common belief that the shocks in the circumstellar medium, pro- 
duced by the violently expanding SNRs, should accelerate the CRs in the Galaxy 
up to about the knee in the all-particle energy spectrum at several 10^^ eV. 

From an experimental point of view it is obvious that any real test can 
only come through characteristic radiation signatures from the dominant nuclear 
particles in individual objects, at comparable energies. The most direct way to 
achieve this test is the detection of gamma-rays at TeV energies. 

2.1. Practical beginning of the source search 

This effort began in earnest about ten years ago when time- dependent nonlinear 
acceleration models in SNRs, still in a hydrodynamic approximation also for the 
CR component, were used to calculate the expected 7-ray emission. The models 
had to assume the form of the energetic particle spectrum, consistent with the 
calculated time evolution of the total energy in CRs, and in detail they turned 
out to be sensitive to the assumed injection rates at suprathermal energies [2,3]. 
While these dynamic models concentrated on the nuclear 7-rays from 7r°-decay 
following inelastic collisions with nuclei from the thermal gas, there was also a 
class of kinematic models by e.g. Naito & Takahara [4] and Gaisser et al. [5] 
who rather assumed distributions of accelerated nuclei and energetic electrons 
but then concentrated on the various radiation processes from these particle pop- 
ulations and their relative importance. Subsequently 7-ray emission models were 
developed from kinetic theory by Berezhko & Volk [6,7]. They involved numerical 
solutions of the full time-dependent, nonlinear CR transport equation in spherical 
symmetry which had been first obtained by Berezhko et al. [8,9]. Based on sta- 
tionary plane wave solutions of the same transport equations with a Monte Carlo 
code [10], Baring et al. [11] also estimated the time- dependent 7-ray emission. 

The conclusions from this theoretical work were cautiously optimistic: nu- 
clear TeV 7-rays from bright nearby objects should be marginally detectable by 
existing ground-based instruments. 

2.2. Early observational attempts 

The first specific 7-ray observations were made in the early nineties. However 
the initial results were rather inconclusive. Esposito et al. [12] had observed sev- 
eral shell-type SNRs, notably G78.2+2.1 (7-Cygni) and IC443, with the EGRET 
instrument on CGRO at 7-ray energies below 1 GeV. Subsequent TeV observa- 
tions by the Whipple [13] and HEGRA [14] groups yielded only upper limits. For 
7-Cygni and IC 443 these upper limits were below the power-law extrapolations 
of the EGRET spectra and, like the Whipple upper limit for Tycho's SNR, they 
tended also to lie below the theoretical estimates [3] for the vr^-decay 7-ray flux. 
In contrast a first TeV-detection was reported for the Southern Hemisphere rem- 
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nant of SN 1006 by the CANGAROO collaboration [15]. The flux value given 
significantly exceeded the expected vr^-decay flux [16]. 

The failure to detect TeV emission from the EGRET sources could be 
understood if in reality a straight extrapolation of the power law spectrum from 
EGRET energies to the TeV range was not required. The expected emissions 
would then come from different objects. Alternatively the remnant could be 
"old" in an evolutionary sense, having lost many of its very high energy particles 
already (see the section on Gas A). Such doubts were substantiated by Brazier et 
al. [17] who rather found a Pulsar in the EGRET 95 % error box for 7-Gygni, 
and Pulsars are generally assumed to exhibit a cutoff in the emission spectrum 
at some tens of GeV. In IG443 Keohane et al. [18] found only two regions with 
hard X-ray emission from high energy (~ 10 TeV) electrons, certainly not over 
most of the radio shell. An association between the hard X-ray sources and the 
EGRET source can not be firmly established either (see [19], which also contains a 
summary of the recent X-ray results). More generally speaking, no unambiguous 
EGRET detection of a shell-type SNR exists. This is one of the questions which 
the GLAST mission will have to resolve. 

2.3. SNR detections in nonthermal hard X-rays and TeV 'y -rays 
In the mid-nineties SN 1006 was detected in hard X-rays [20,21] and much of 
this emission was attributed to synchrotron radiation of electrons with energies 
of tens of TeV [20,22], see Fig. 3a. Subsequently also several other shell-type 
remnants where shown to have power law tails in their emission, notably SNR 
RX J1713.7-3946 in the Southern Hemisphere [24,25] and Gas A [26]. 

The question was then, whether the energetic electrons would not also be 
visible in TeV 7-rays through their Inverse Gompton (IG) scattering on ambient 
low energy photons, primarily from the Gosmic Microwave Background (GMB) 
[27,28]. In 1998 the GANGAROO collaboration published its detection of SN 1006 
[15] mentioned earlier, and interpreted it in terms of IG emission. A subsequent 
TeV detection of SNR RX J1713. 7-3946 with the same telescope [29] was given an 
analogous interpretation. The only critical discussion at the time was published 
by Atoyan and Aharonian [30]. 

The prevalent opinion that there was no need for nuclear GRs to explain 
the TeV 7-ray emission from SNRs gave room to doubts as to the existence of 
any significant quantities of GR nuclei in shell- type SNRs at all [31]. However this 
pessimistic turn was provoked by phenomenological arguments not by theory, as 
all the arguments before. Again on a phenomenological basis, the GANGAROO 
collaboration on the other hand reversed its view regarding SNR RX J1713.7- 
3946 by presenting new arguments which now favored a hadronic 7-ray origin [32] 
(see also [25] above). This led to a controversial discussion, where several groups 
[33,34] gave empirical counterarguments which questioned the new interpretation. 
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Fig. 3. a. Synchrotron and 7-ray emission from SN 1006. Radio data are from [22], 
X-ray data are from [21] and [23]. Approximately a power law plus cutoff (dashed 
curve) for the electron distribution in a weak field of several iiG has been used in 
a phenomenological fit of the data. b. The same electron spectra are used for an 
IC fit to the differential 7-ray energy flux data reported in [36]. Also the EGRET 
upper limits are shown. 



Unfortunately SNR RX J1713. 7-3946 is a complex source which may either be the 
result of the thermonuclear explosion/deflagration of an accreting White Dwarf 
(SN Type la) that ejects a Chandrasekhar mass, or of the core collapse of a 
massive star. The distance uncertainty ranges from 1 to 6 kpc [24, 25]. The 
remnant is situated in a complex interstellar environment (e.g. [35]) and is far 
less well studied than SN 1006, for example. Therefore the debate remains quite 
inconclusive at present. SN 1006 is a much clearer case. It is a SN type la, with 
excellent radio and X-ray observations and extensive morphological studies. Fig. 
3b shows the latest published 7-ray results [36], with an IC fit for the spectrum. 

Our own conclusion was that phenomenological considerations are impor- 
tant but that they will not give a convincing answer to the question of the accel- 
eration of CR nuclei in SNRs as reasonable as they may appear individually, one 
by one. There is a need for quantitative comparison of the observations with a 
consistent theory, as in other areas of physics. Contradictory interpretations of 
the observations can be excluded only if such a detailed picture is available. 

3. Theory of diffusive shock acceleration as applied to SNRs 

Let me summarize such a theory in three cartoons. Obviously they are only meant 
to illustrate the essential physics; the real description is given by the equations 
which I shall only write down to hint at the mathematical aspects. Much of the 
theory is reviewed in e.g. [37,38,39,40,41]. 

The left cartoon of Fig. 4 shows the spatial dependence in a plane shock, propa- 



7 



In f 




Fig. 4. Spatial dependence of flow velocity Ui, gas pressure (denoted here by p), 
mass density p, and CR pressure pc, in the regions upstream (i=l) and downstream 
(i=2) of a plane shock, as a function of the space coordinate x in the shock frame of 
reference. The wavy line indicates the growing amplitude of the scattering magnetic 
fluctuations of the Alfvenic wave field which is excited by the diffusive CR current 
—ndpc/dx into the upstream precursor region and compressed in the subshock at 
x = 0. Vi is the wave phase velocity vector in the upstream fluid frame. Particles 
gain energy by repeated scatterings across the shock (left panel). The right panel 
shows the momentum dependence of the downstream particle distribution. The 
three regions joined by dashed lines correspond to the thermal plasma, characterized 
by the thermal energy (-Eth) and extrapolated by a dotted curve, the low-energy 
nonthermal particles (Eiow), and the high-energy particles (Ehigh), respectively. 



gating parallel to the magnetic field. The nonlinear backreaction of the accelerated 
particles on the dynamics of the thermal gas through their pressure gradient cre- 
ates a smooth precursor, followed by a sharp subshock, where the fluid quantities 
of the thermal gas jump disco ntinuously. The CR pressure 

Pc = ^ dp / 2, 2 f 

O J Pin J VP + ^ C 

is given in terms of the isotropic part of the particle distribution function /(x, p, t), 
averaged over the magnetic fluctuations; the quantities p and m denote the mag- 
nitude of particle momentum and its mass. In the downstream distribution (right 
panel) CRs are characterized by momenta above the region of injection where 
downstream particles can outrun the shock along the magnetic field. The non- 
thermal spectrum starts with a steep power law distribution, characterized by 
diffusive acceleration at the subshock only. High energy particles " see" the entire 
compression as a discontinuity and have therefore a harder spectrum that begins 
in the region where protons become ultrarelativistic [9]. 

Strong shocks with Alfvenic Mach numbers Ma = Ui/Vi ^ 1 typically 
inject so many suprathermal nuclear particles into the acceleration process that 
the acceleration is efficient, Pc ~ PiUf, and that nonlinear backreaction leads to a 



8 



strong precursor. In SNRs Ma is of order 10^ at very early times, and for SN 1006 
it is still as large as 150 at the present epoch. Thus the total compression ratio 
Ttot > 4, whereas the subshock compression ratio rgub < 4. Only part of the 
overall energy dissipation goes into thermal energy at the subshock, the rest is 
given to the energetic particles and to the growing waves. 

The most difficult aspect of particle acceleration is the magnetic field and 
its fiuctuations, because this field is what the energetic particles interact with 
directly. And the field is not simply given externally, but is rather excited by the 
particles themselves as another part of the nonlinearity of the process. 

The selfconsistent wave production by the accelerated particles themselves 
can only be estimated rather approximately. A simple-minded quasilinear calcula- 
tion [42] for efficient acceleration leads to an extremely large wave field 6B/B ^ 1, 
where 6B is the total rms wave amplitude and B is the assumed mean magnetic 
field strength. This conclusion has been shown to hold also differentially for the 
wave energy per unit logarithmic bandwidth in wave number in relation to the 
resonant energetic particle pressure Pc{p) = 47rc/3 f{x,p,t) / y/p^^^rwFc? per 
logarithmic interval in p [41]. 

However, such a result defeats the original assumption SB/B ^ 1 of per- 
turbation theory, on which the derivation of the result is predicated (I will loosely 
continue to use the term 6B for both the differential and the integral wave am- 
plitude). The question is then how to deal with this outrageous situation. In 
fact, theory has no precise solution at the moment. Apart from wave damping, 
the simplest and physically most plausible conclusion is that the strong wave 
production amplifies the mean field to an effective field B^ff ^ Bq [38], with 
6B/Beff ~ 1. In this limit the scattering mean free path Xmfp reaches its mini- 
mum value Tgyroip), where rgyro is the particle gyro radius in the field -Be//- This 
is the so-called Bohm limit in the effective field and implies an isotropic diffu- 
sion particle coefficient i?c//) = vrg{p)/3. Field amplification has recently 
been successfully calculated in a simplified nonlinear model of wave turbulence 
by Lucek & Bell [43] and Bell & Lucek [44]. The result is the Bohm limit. As a 
scaling relation for strong shocks it is already suggested by perturbation theory. 

Finally strong wave production has also profound consequences for the rate 
and the geometry of ion injection (third cartoon of Fig. 5a): where on the shock 
surface injection can occur in the first place - in an approximately spherical shock 
propagating in an inhomogeneous external magnetic field this is only possible in 
the "polar" regions - it occurs only at those field lines that are instantaneously 
quasi-parallel to the subshock normal direction. On all other field lines injection 
is suppressed until the field changes back to quasi-parallel. This has the conse- 
quence that the ion injection rate (the fraction of the incoming fiux of thermal 
particles that can be accelerated) is reduced by about two orders of magnitude 
compared to that in a laminar magnetic field. The precise numerical value of the 
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reduction factor depends on the wave spectrum. For a uniform ambient field the 
systematic reduction of injection will lead to dipolar asymetries in the nonthermal 
morphology and requires a renormalization of the overall particle production rate, 
as calculated in spherical symmetry, by a factor fre < 1 [45]. 

The resulting theory, applied to the dynamics of a SNR and assuming 
spherical symmetry, looks as follows: 

CR transport and acceleration at shocks is described by a Fokker-Planck 
equation for the isotropic part of the mean particle distribution fi{p,r,t) for ions 
(p) and electrons (e), respectively: 

— - V/iV/i + f/V/i + — /i = <5i' (2) 

ot 6 op p"^ Op \Ti J 

where the coefficients in the transport eqs. are given by the mean mass velocity 
U{r,t) and its divergence, as determined from the hydrodynamic eqs., and by a 
common diffusion coefficient k{p, B^ff) as a result of scattering by the fluctuating 
magnetic fleld. Qi denotes the injection rate, and Ti is the loss time, taken to be 
zero for ions and equal to the synchrotron loss time for electrons. 

The hydrodynamics of the thermal plasma, even if we consider it as an 
ideal fluid, is in turn nonlinearly coupled with the energetic particles through the 
CR pressure gradient dpc/dx, deflned above, as well as wave dissipation, see e.g. 
[7]: 

^ + V{pU) = 0, (3) 

dU - - 
p— + p(UV)U = -Vip, + Pg), (4) 

^ + iUV)pg + lgiyU)pg = - 7,)caVpe, (5) 

where p, 7^ = 5/3 and Pg denote the mass density, gas speciflc heat ratio and 
gas pressure, respectively. Gas heating due to Alfven wave dissipation in the 
upstream region is described by the parameter < 1. 

Already in this lowest approximation the theory involves a time-dependent 
nonlinear system of coupled partial integro-differential equations. 

In principle the theory is quite general and can be used to describe any 
ionized system in which matter dominates and thermal gas motions are nonrel- 
ativistic. The above application to point explosions like SNRs is the simplest 
one from an analytical point of view. Astrophysically speaking it is nevertheless 
fundamental. In this case the theory contains only two "parameters" that are not 
very well calculable quantitatively: -Be//? and the ion injection rate Qp. They are 
strongly related through the nonlinear development of wave production discussed 
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t/to 

Fig. 5. Left panel: Ion injection into diffusive shock acceleration. Only for @nB2 suffi- 
ciently smaller than 90° suprathermal ions can escape upstream along the magnetic 
field 5. In a spherical SNR in a uniform field B^i, injection, strong wave production 
and acceleration is only possible in the hatched polar regions. Therefore hadronic 
7-ray emission is dipolar and the same is true for the synchrotron emission as a 
result of field amplification. The right panel contains the temporal evolution of the 
overall dynamics for SN 1006 [46]. For external density n = 0.3 cm~^ the shock 
radius Rs and velocity Vs are exhibited together with the observed values. Total, a, 
and subshock, ag, compression ratios are given in the middle panel. The quantities 
EsN, Ef>j, Egk, Egt and Ec denote the total SNR mechanical, ejected, gas kinetic, 
gas thermal, and un-renormalized CR energy contents, respectively. 
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above. It is therefore plausible to use an experimental/observational input for 
these parameters in order to be free from uncertain theoretical approximations. 

How do we achieve this observational input for SNRs? The natural solution 
is given in terms of the electron component and the corresponding sychrotron 
observations [46]. The electrons are parasitically accelerated together with the 
nuclear particles which generate the scattering wave field through their dominant 
mass and energy density, at least above the (unknown) electron injection energy. 
This means that above injection the overall electron momentum distribution must 
have equal form to that of the nuclear particles except for radiative losses. The 
synchrotron observations then yield three quantities: (i) the radio synchrotron 
spectrum will be steeper than in the test particle approximation; interpreting this 
in terms of nonlinear modification determines then the ion injection rate (ii) the 
requirement that these radio electrons have energies < 1 GeV gives the value of 
-Be// (iii) the electron density amplitude from the magnitude of the synchrotron 
emission determines the electron:proton ratio. The actual average nuclear CR 
density in the SNR is reduced from the calculated spherically symmetric ion 
amplitude by the ratio fre of the area of efficient injection to the total SNR 
surface area (Renormalization) [45] . 

4. SN 1006 

In the radio and in nonthermal X-rays SN 1006 has a dipolar structure [20,47]. 
This is consistent with ion injection theory, as described above and with a diam- 
eter of 0.5° it should be an extended TeV source as well, allowing a test of the 
7-ray morphology. From astronomical measurements at radio and X-ray wave- 
lengths the distance is about 1.8 kpc, and the expansion velocity amounts to 3000 
km/sec. A rather uncertain astronomical parameter is the density n of the am- 
bient Interstellar Medium (ISM). In the following we use a value n = 0.3 cm~^ 
which is suggested by X-ray measurements. However, the density could be as low 
as n = 0.1 cm~^ since SN 1006 is located relatively far above the Galactic Plane. 
This also suggests a standard interstellar magnetic field of several /iG, rather uni- 
form on the scale of the SNR. However, as emphasized below, the effective field 
is expected to be substantially higher than this. 

4.I. Model calculations for SN 1006 

I shall now paraphrase the detailed model calculations [46] that are based on the 
theory discussed before. The overall hydrodynamic quantities are given in Fig. 
5b. From a simultaneous fit of the observed SNR radius and expansion speed at 
the known age of the system of almost 1000 years the total mechanical explosion 
energy must be chosen as Esn = 3 10^^; Esn is somewhat high compared to the 
canonical 10^^ erg. The total shock compression ratio is substantially above the 
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adiabatic value for a strong shock. The spherically symmetric calculation gives a 
total CR energy E^. !^ O.GEsn at late times. However, the renormalization factor 
is 0.2 [45]. This results in a fraction of about 10 % in nuclear CR energy at the 
late Sedov phase, when the source particles are expected to be released into the 
ISM. 




Fig. 6. Volume-integrated particle distributions N and synchrotron spectra I'Si, for 
accelerated protons (p) and electrons (e) in SN 1006 [46]. The observational data 
are the same as in Fig. 3a. 

The spatially integrated spectra of CR electrons and nuclei (Fig. 6a) harden to- 
wards the cutoff, except for electron cooling. The high effective internal field 
-Be// ~ 120 /iG is determined from the spectral slope of the radio synchrotron 
spectrum in Fig. 6b which also fixes the injection rate to about 10~^ times the flux 
of incoming particles. The flattening of the high frequency synchrotron spectrum 
also allows a smooth connection below the thermal X-rays to the hard X-ray data 
points. All this implies a somewhat low predicted electron:proton ratio in the 
overall number of accelerated particles of about 2 ■ 10~^ for this source, compared 
to the canonical interstellar value of 10~^. 

The predicted contributions to the differential 7-ray energy spectrum, inte- 
grated over the SNR volume, are shown in Fig. 7a and the radial 7-ray brightness 
at 3 TeV is exhibited in Fig. 7b. For the assumed ISM density n the IC con- 
tribution from CR electrons is about one order of magnitude below the hadronic 
contribution due to Tr^-decay, which in turn is still below the EGRET upper limits. 
As consequence of synchrotron cooling, the IC emission spectrum has an almost 
identical form to the hadronic 7-ray spectrum from several GeV up to the cutoff. 
This makes attempts for a distiction of spectral slopes very difficult except at 
very low and very high energies. 

It is in this context important to note that the vr^-decay 7-ray flux scales 
like ra^ to lowest order, although the strong nonlinearity of the system makes this 
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scaling only a rough first approximation. Decreasing the external density by a 
factor of 2 will therefore diminish the hadronic 7-ray flux by a factor of about 4. 
At the same time the electron:proton ratio will increase, whereas the derived total 
energy Esn will decrease. An external density which is lower than n = 0.3 cm~^ 
may explain the fact that the individual H.E.S.S. telescopes in Namibia have not 
been able to detect SN 1006 at the level suggested by the CANGAROO results. 

Fig. 3b in fact indicates an IC spectrum similar to the one used to explain 
the more recent data points given by the CANGAROO II telescope [36]. Such 
a spectrum is the result of a scenario that neglects acceleration of nuclear par- 
ticles and assumes an average internal field of about 10 /iG, appropriate for an 
MHD-compressed ISM magnetic field of several fiG. The form of this spectrum is 
different from both the hadronic and the cooled IC spectra above a few GeV. Nev- 
ertheless, a well-defined observational 7-ray spectrum is required to distinguish 
this scenario from the solid line spectra calculated from the nonlinear theory if 
the fluxes are comparable. 




Fig. 7. Predicted differential energy spectrum of vr'^-decay and IC emission (left 
panel), and 7-ray brightness (at 3 TeV) as functions of radius r in units of the 
shock radius Rs (right panel), for SN 1006 for an ambient density n = 0.3 cm~^. 
Upper limits from EGRET and the reported CANGAROO fluxes are indicated [46]. 



The picture becomes richer when we consider the morphology as well, for 
instance the radial dependence in the neighborhood of the dipole axis (Fig. 7b, for 
a 7-ray energy of 3 TeV). Here the IC emission is concentrated into an extremely 
small scale le ~ 10~^Rs near the shock at r = i?s as a consequence of synchrotron 
cooling of the generating ~ 100 TeV electrons. Also the hadronic component is 
fairly strongly confined to a narrow shell, since the thermal gas is concentrated 
there. The emission morphology for a low-field scenario (not shown here) would be 
considerably more extended, only exhibiting the large-scale adiabatic expansion 
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losses in the interior. In summary, the 7-ray emission is expected to show a 
dipolar spatial morphology [48], similar to that of the synchrotron emission. 

An important observational result has been recently been found in Chandra 
observations at X-ray energies in the 1 to 10 keV region [47,49] (Fig. 8a). 




Fig. 8. Northeast limb of SN 1006, as seen by Chandra between 0.5 and 2 keV [47] 
(left), and projected X-ray brightness of the most narrow structure, compared with 
the theoretical prediction [50] (right). Vertical dashes are from [49], histogram from 
[47]. The thin curves correspond to the low-field scenario in Fig. 3a. 

They demonstrate very narrow spatial X-ray emission structures, on a scale of 
10~^-Rs; confirming the magnetic field amplification made possible by a nonther- 
mal particle population dominated by nuclear particles [46,50]. In contrast, the 
low-field scenario clearly fails when compared to the observations (Fig. 8b). 

As a corollary, the electron population cannot create the strong field fluc- 
tuations and the required high effective field by itself [50]. This means that the 
dominant accelerating nuclear component is not only possible experimentally as 
well as theoretically: it is also necessary. 

5. Cassiopeia A 

This SN type lb, whose progenitor was probably a massive Wolf-Rayet star, is 
on the other end of the scale of Supernovae when compared with SN 1006. Such 
massive progenitor stars modify the circumstellar medium through substantial 
mass-loss during several successive stehar wind phases: the fast rarefied Blue 
Supergiant wind during the main sequence phase turns into a massive but slow 
Red Supergiant (RSG) wind at a late stage, to be subsequently compressed from 
inside by yet another fast wind, now from the emerging Wolf-Rayet star, until this 
star finally collapses as Supernova. The turbulent compressed shell of the RSG 
material has a high gas density of about 10 cm~^ owing to radiative cooling and 
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is identified with the so-called bright ring in Cas A [51], see Fig. 2. According to 
this picture, the SNR's leading shock has already reached the unperturbed RSG 
wind region. 




Fig. 9. Schematic of Cas A's circumstellar environment at different stages of evo- 
lution. Left panel: At explosion the turbulent shell is bordered by the interior 
Wolf-Rayet wind bubble and the exterior Red Supergiant wind region. Right panel: 
Today's SNR shock and clumpy ejecta (white spots) beyond the shocked shell and 
part of the Red Supergiant wind region, superposed on radio image of Fig. 2. 
(Courtesy G. Piihlhofer). 

The HEGRA telescope system observed Cas A for 232 hours, the longest 
pointing used in 7-ray astronomy until now, and finally detected the source at 
3.3 percent of the Crab level [52]. The acceleration model summarized here [53] 
basically follows the picture described above and assumes a Parker spiral type 
mean circumstellar magnetic field topology. The initial configuration before ex- 
plosion is schematically shown in Fig. 9a, while Fig. 9b pictures the situation of 
today. 

Given the synchrotron spectrum, the evolution of this complex system 
containing the turbulent wind shell allows a consistent description of the SNR 
dynamics with a total mechanical energy Esn = 4 ■ 10^° erg, a present shock 
speed Vs = 2000 km/sec, and a very large downstream effective field -Be// = 1 niG 
in the shell. Even though the available synchtrotron observations are seperated 
by 30 years, about 10 % of the remnant's age, a good fit to the hardening radio 
spectrum is possible, followed by a synchrotron cooling range in the infrared to 
X-ray range, ending in a cutoff at nonthermal hard X-rays (Fig. 10a). Also the 
observed temporal decrease of the synchrotron emission is consistent with the 
picture of shock propagation into an expanding wind region of decreasing gas 
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Fig. 10. Calculated synchrotron spectra at epochs 1970 {solid curve), 2002 (dashed 
curve) and 2022 (dashed- dotted curve), respectively, together with observational data 
over a comparable time span (left panel). The right panel shows the predicted inte- 
gral ^ -ray energy spectrum: Nonthermal Bremsstrahlung (NB, dashed line) and IC 
emission (dash-dotted line), together with dominant -k^ -decay emission (solid line) 
are far below the EGRET upper limit, but agree with the HEGRA flux at 1 Te V. The 
dotted line indicates a possibly lowered proton cutoff as a result of particle escape 
from the SNR [53]. 

density. 

The predicted differential 7-ray energy spectrum is shown in Fig. 10b. 
As a consequence of the high values of gas density and effective magnetic field 
strength, it is completely dominated by hadronic 7-rays from vr^-decay which are 
still two orders of magnitude below the EGRET upper limit. The hadronic flux 
has been renormalized by a factor f^e = 1/6 [45]. This prediction reasonably 
agrees with the observed 7-ray flux value around 1 TeV. The predicted IC and 
nonthermal Bremsstrahlung fluxes are very low. We therefore conclude that Cas A 
is a hadronic 7-ray source. The predicted proton spectrum reaches the knee 
region. Even if escape of CR nuclei from the weakening shock is likely to set in 
already at this low age, Cas A could be considered as a member of the hypothetical 
CR source population of Galactic SNRs. 

6. A few personal theses 

1 want to summarize this talk in a few personal theses. I call them personal 
because they may not reflect the position of the two other speakers. On the 
other hand, the conclusions are based on the joint efforts of E.G. Berezhko, L.T. 
Ksenofontov, G. Piihlhofer and myself, and to this extent they are not only my 
own conclusions: 

• The astrophysical processes of CR origin up to the knee are basically un- 
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derstood. The sources are the Galactic SNRs. This concerns 99.9 percent 
of the total energy density in CRs. 

• Gamma-ray observations for SNRs are critical since only 7-rays have in- 
dividual energies comparable to those of the generating charged particles. 
Currently experimental results are still scarce. 

• The nonthermal synchrotron properties are an indispensible multi-wavelength 
aspect of SNRs. This concerns both the morphology and the spectrum. 

• On account of the hard SNR spectra, the most appropriate instruments 
at 7-ray energies are large ground-based detectors like CANGAROO III, 
H.E.S.S., MAGIC and VERITAS. They will play a decisive role for the de- 
termination of the morphologies as well as the spectra of individual objects. 
They are needed to increase the detection statistics of SNRs. 

• The most important global astronomical test is the detection and spectral 
decomposition of the " diffuse" TeV 7-ray background in the Galactic Plane 
to which the entire Galactic SNR population contributes. 
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